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Abstract. We thoroughly investigate vibrational quantum dynamics of dimers 
attached to He droplets motivated by recent measurements with K2 [T]. For those 
femtosecond pump-probe experiments, crucial observed features are not reproduced by 
gas phase calculations but agreement is found using a description based on dissipative 
quantum dynamics, as briefly shown in P]. Here we present a detailed study of 
the influence of possible effects induced by the droplet. The helium droplet causes 
electronic decoherence, shifts of potential surfaces, and relaxation of wave packets in 
attached dimers. Moreover, a realistic description of (stochastic) desorption of dimers 
off the droplet needs to be taken into account. Step by step we include and study 
the importance of these effects in our full quantum calculation. This allows us to 
reproduce and explain all major experimental findings. We find that desorption is fast 
and occurs already within 2 — 10 ps after electronic excitation. A further finding is 
that slow vibrational motion in the ground state can be considered frictionless. 
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1. Introduction 

Helium nanodroplet isolation (HENDI) spectroscopy allows to study atoms, molecules 
and clusters embedded in an ideal cryogenic environment [3] . Ultracold helium droplets 
(0.38 K) provide a gentle, since weakly disturbing host for embedded species, which 
can be studied with high resolution spectroscopy [1]. However, spectra of immersed 
species are slightly broadened and shifted away from their gas phase value due to 
their interaction with the surrounding He droplet. Inhomogeneous broadening in pure 
rotational spectra is ascribed to the motion of the purity inside the droplet O |6] or 
to coupling to collective degrees of freedom of the droplet [7]. The overall spectral 
features, which are unseen in gas phase measurements, have been used to study the 
weak interaction between dopant and helium droplet [HI El ttHl El [HI [12] • Attached 
species allow to probe the peculiar properties of the superfluid He droplet itself [13]. 
Further insight into the interaction and the quantum properties of the droplet is gained 
from recently obtained time- resolved studies [H [21 [H [IS] • 

Helium nanodroplets, typically consisting of several thousand ^He atoms, are ideally 
suited to study relaxation (cooling) of embedded species [IHl HZl UHl UHl ISO]- Whether 
dissipation plays a role depends on the involved energy scales, the couphng strengths, 
or on typical time scales in the system and "bath" [211 [22] . The group of Miller has 
studied relaxation of the HF system inside He droplets [23] . They report on ineffective 
vibrational relaxation due to a mismatch in energy scales. However, rotational relaxation 
in the immersed HF system is fast and appears as Lorentzian line broadening in the 
rotational transition spectrum, (see also [7] and references therein). 

Alkali metal atoms and molecules are known to reside in bubble-like structures on 
the surface of He droplets [211 [251 [26] . Attached Rb2 dimers in the triplet state reveal the 
presence of vibrational relaxation on the measurement timescale [23 US]. For Lithium 
dimers interacting with a He environment, vibrational relaxation has been investigated 
by means of Monte Carlo calculations p8j, including a few He atoms. The corresponding 
relaxation rates depend on whether the dimer is orientated in-plane or out-of-plane with 
respect to the He surface [29] . 

Femtosecond pump-probe techniques are established tools to analyze the ultrafast 
vibrational motion in molecular systems [301 IHH [32] • A first laser pulse excites a coherent 
wave packet, which is allowed to freely evolve on the respective energy surface. The WP 
is probed by a time-delayed pulse. Since the early studies by the Zewail group [531 [M], 
vibrational wave packets in various diatomic systems have been studied, such as Na2 
[53] or K2 [Sni [571 [3H] • When molecules are located in a solid rare-gas matrix [511 [10] , a 
suppression of revival structures in the pump-probe signal indicate a loss of vibrational 
coherence. Meier et. al. have thoroughly investigated such decoherence in molecular 
systems placed in a rare-gas environment [H], motivated by seminal experiments in 
the Zewail group [12]. For experiments at room temperature with a significant thermal 
occupation of rotational states, unavoidable coupling between internal vibrational and 
rotational degrees of freedom may lead to sufficient decoherence to suppress (fractional) 
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revivals [13] . Experiments with Rb2 dimers on He droplets reveal an ongoing decay of the 
pump-probe ion yield |T1] . The decay has been ascribed to damping and accompanying 
decoherence of vibrational wave packets [15]. In this work we investigate potassium 
dimers on He droplets, studied experimentally with the pump-probe technique [1]. A 
brief account of a theoretical description based on dissipative quantum dynamics was 
given in [2]. Experimental spectra show significant deviations from corresponding gas 
phase calculations for the unperturbed dimer. 

In this paper, in a phenomenological approach, we investigate in detail how the 
helium influence may be described effectively. First, we see that the helium environment 
destroys electronic coherence, which is imprinted by the exciting laser pulse. As possible 
causes for electronic decoherence, we consider a distribution of shifts of electronic 
surfaces. Electronic decoherence alone cannot account for the decay of the signal 
at certain excitation wavelengths. Therefore, in a next step, we include a general 
damping of vibrational wave packets on each electronic surface. The effective dynamics 
is described by means of a quantum optical master equation. 

Moreover, it is important to take into account the desorption of dimers from the 
droplet properly. No general rule can be given when attached atoms or molecules 
leave the droplet. It is known that lighter alkali metal atoms leave the droplet upon 
electronic excitation [H] or form a bound exciplex [33]. The exciplex tends to desorb 
off the droplet surface during the formation process or several picosecond thereafter 
[46j. On the other hand, desorption of Rb atoms may be completely inhibited upon 
electronical excitation in a certain laser wavelength range [1?]. Recent measurements 
indicate that K2 molecules desorb several picoseconds after laser excitation [1] . Indeed, 
(stochastic) desorption of dimers are a crucial ingredient to explain spectral features 
with a theoretical model [2] . 

We here report on a more realistic description of desorption - we use a model which 
only allows electronically excited molecules to leave the droplet. We thus extend our 
previous, state- independent desorption scheme [2]. Together with dissipation and shifts 
of potential energy surfaces, one can explain experimental findings over the full laser 
excitation range reported in [1]. 

The article is organized as follows: In section [2] we review the calculation of the 
pump- probe signal for free dimers. In the following section [3} the influence of the helium 



environment is considered in a phenomenological way. In section 3.3 we explain how we 
treat desorption of dimers. Finally, for an even better agreement with experiment, we 
consider undamped motion in the electronic ground state. A thorough comparison 
between theoretical and experimental findings is given at every step in the model. 
Section |4] is devoted to the conclusions. 

2. Pump-probe signal 

A first "pump" laser pulse creates a vibrational wave packet (WP) {ipi) in some electronic 
state \i). The excited WP oscillates in a region between classical inner and outer turning 
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Figure 1. Excitation schemes in the potassium dimer for two distinct laser 
wavelengths A — 833 nm (scheme I) and A — 800 nm (scheme II). While the first 
scheme exclusively maps the WP in the A^E+ state, the second scheme allows to 
observe the WP in the ground state X-'^S+. 



point of that surface. It may periodically enter and leave a transition region, where a 
resonance condition with higher lying states is met and the potential energy difference 
matches the second, "probe" pulse energy. This region defines the so-called Franck- 
Condon (FC) window. The time-delayed probe pulse leads to a significant number 
of ions only when the WP is located in the FC window. By varying the time delay r 
between pump-and probe pulse, one obtains an oscillatory ion yield. A typical excitation 
scheme is depicted in fig. [l](a), where one probes the dynamics of the WP on the A ^Sj^ 
surface. 

Wave packets may be excited in several involved electronic surfaces of the dimer. 
More specifically, the pump pulse prepares the dimer in a superposition of electronic 
states, such that the full state vector |^) takes the form 

|^(W)) = $^|^.)K). (1) 

i 

Here, ipi denotes the WP on a specific electronic state i and 

P^={i^^\i^^) (2) 

is the probability that the electronic state i is excited. 

For the pump-probe signal, we fully numerically solve the time-dependent 
Schrodinger equation (TDSE) 
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for the full state vector (see also |1H 



)]). In the matrix equation, the diagonal 
elements denote the molecular Hamiltonian -ffmoi = Ylii Hi = T + J2i ^5 which involves 
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kinetic energy T = (reduced mass /i) and adiabatic potentials Vi. In the final 

ionic state, the energy E of the ejected electron is included in the diagonal entry of the 
Hamiltonian H^^e^ such that Vi=finai = V3+E. The coupling to the laser field is described 
by the matrix elements Jij = —fIij-E{t), where /lij denotes the transition dipole moment. 
Both pump- and probe pulse have the form -Epump/probe(^) = eo£(t) cos {oJi,t) and ujl is the 
respective laser frequency. Moreover, eo is the polarization and e{t) the shape function 
of the field, which is assumed to be Gaussian. For the field parameters we use a full 
width at half maximum of the laser pulse of llOfs and an intensity of 1.2 GW/cm^. The 
employed intensity is higher than the experimentally estimated value (1=0.5 GW/cm^ 
PP), but, according to still located in the moderate power regime. The ionic state 
ip3{E) also depends on the energy E of the ejected electron. We use discretization of the 
(electronic) continuum, a technique successfully employed earlier [IH]- One determines 
the final state probability \ip3{Ek)\^ after both pulses have passed for distinct electronic 
energies E^. The pump-probe signal is proportional to the sum over different electronic 
contributions. 



In the limit t — )■ 00 (upon complete decay of the second pulse), it only depends on 
the delay r between the pulses. For the propagation of the wave function, we use the 
split-operator method [lU] . 

The ion signal is composed of beat frequencies u^v' = {E^ — Eyi)/h between all 
pairs of energy levels that contribute to the WP [T5,|. The most prominent oscillation 
originates from the energy spacing between central and neighboring vibrational levels v 
and V ±1. This oscillation has the frequency 



and is characteristic for the electronic surface i. In an anharmonic potential, the 
level spacing and therefore also cuj decreases as v increases. Higher-order frequency 
components Wi},t}+At; with Af > 1 are visible in the Fourier transform (FT) of the signal. 

The laser wavelength A determines which wave packets ipi can be mapped to the 
final state. In the one-color pump-probe setup, we consider two different excitation 
schemes. 

For 820 nm ^ A < 840 nm (scheme I), one exclusively follows the vibrational 
dynamics in the excited state A^S+. The WP in that state can be probed at the 
outer turning point (see fig. [l]^a)). Contributions from other surfaces are negligible. 
For this scheme we concentrate on an excitation wavelength A = 833 nm. Fig. [2]^a) 
shows the theoretical gas phase calculation. Fig. [2]^b) shows experimental spectra 
at this wavelength obtained from dimers attached to helium nanodroplets pp. In 
the experiment, the signal amplitude significantly decreases, but oscillates with nearly 
constant amplitude at later delay times. 

For 800 ^ A < 820 nm (scheme II), transitions preferably take place at the inner 
turning point of the excited WP. The wave packet in the electronic ground state X 
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Figure 2. Pump-probe ion yield at A = 833 nm. (a) Numerical gas phase calculation, 
(b) Experimental result (from [T]). The oscillation is exclusively attributed to the 
circulating WP in the A-^E+ state. 



is excited through resonant impulsive stimulated Raman scattering (RISRS). It can be 
probed through a 3-photon process, see fig. [1} On the other hand, simultaneous and 
coherent contributions from WPs in electronic excited states lead to constructive or 
destructive interferences and thus to an unstructured ion yield pS]- We use sliding 
window Fourier transforms J^{u}, r) to follow the evolution of respective beat frequency 
components as a function of the delay time [501 EI] • Fig- shows a window transform 
(spectrogram) for excitation at A = 800 nm. A frequency component ua ~ 63 cm~^ 
can be ascribed to the WP in the state A^S+. Moreover, a higher-order frequency 
component around ~ 2uja is visible as is a significant contribution at ~ 85 cm^^ after 
about 15 ps that can be traced back to originate from the 2^Ilg surface. A contribution 
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Figure 3. Spectra 7^(10, t) in the time-frequency domain at A — 800 nm. (a) 
Numerical gas phase calculation, where electronic interferences lead to a spectrum 
that is difficult to interpret, (b) Experimental HENDI result (from 1 ). Electronic 
interferences are absent, allowing for a clear identification of the structures of the 
spectrum. 
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Figure 4. Shown is a the spectrogram r), which is obtained upon projecting 

the created fuh wave function '^{t) on the electronic ground state [i — 0) and (b) 
first excited state (i = 1). (c) In the incoherent sum, interferences are excluded. In 
particular, the ground state component is clearly resolved. 



from the ground state and resulting frequency component ojx^ however, is missing. 
In agreement with earher findings [3H1 136] , the mapping of the ground state WP is not 
possible at low to intermediate intensities in the gas phase. This is because, as explained 
below, contributions from different potential energy surfaces interfere destructively. In 
contrast, for dimers attached to He droplets, the contribution from the state A^E+ is 
suppressed after about 8 ps, while on the other hand the vibrational ground state WP 
is clearly resolved, see fig. |3](b). 

Qualitatively, the difference between fig. [sj^a) (gas phase) and fig. |3](b) (experiment) 
can be explained by the loss of electronic coherence alone: In fig. [4], we construct an 
artificial signal from the incoherent sum of contributions. The full coherent wave vector 
after decay of the pump pulse at t = tpump, eq. ([T]), serves as a starting point. In order 
to determine the contribution of a single WP ipi, we project the fully coherent wave 
vector according to \^i(tpump) = -Pi|^(''^pump)) with the projector Pj = \i){i\ on one of the 
electronic states. The electronic occupation pi is not altered and the usual probe scheme 
is employed after projection. The resulting spectrogram F^u, r) after projection on the 
state A^S^ and X^S^, respectively, is shown in fig. |4](a)/(b). The incoherent sum of 
contributions is given by 

^'-(u;,r) = ^J-^(u;,r). (6) 

i 

Contrary to the coherent signal, fig. |3|^a), in eq. ([6]) interferences are removed by hand 
and contributions from ground- and excited state WPs are clearly visible. Obviously, 
the incoherent sum fig. |4]^c) already reproduces important features of the experimental 
data. In the experiment (fig. |3](b)), however, the component ua fades out at later delay 
times, while the component ux becomes dominant. Therefore, a model based on pure 
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loss of electronic coherence alone cannot explain the measured spectrum. One has to 
consider additional influences of the He environment on the dimer dynamics. 

3. He influence and results 

In a phenomenological approach, we take into account three possible effects: 

(i) A He-induced energetic shift of electronic potential energy surfaces, with possibly 
small fluctuations. 

(ii) Damping of vibrational wave packets. This effect is treated fully quantum 
mechanically within the master equation approach. The WP ipi dissipates energy 
with a certain damping rate 7^. This rate is here seen as a fit parameter and is 
adjusted to experimental observations. 

(iii) Desorption of dimers off the droplet. After desorption, the influence of the helium 
droplet (shift /damping) vanishes. 

A description based on damping of vibrational wave packets has been applied to HENDI 
studies with spin triplet Rb2 dimers [TS]. There, slow vibrational decoherence (as a 
consequence of very weak dissipation) is most relevant for the decay of the revival 
amplitude. In particular, as witnessed by the ongoing decay, desorption of dimers seems 
to be inhibited. 

In the following, step by step, we include shift, damping and finally desorption in 
the calculation of pump-probe spectra. It is found that all model "ingredients" 1.-3. are 
crucial to find agreement with obtained experimental spectra. 

3.1. Electronic shifts 

The He environment may lead to shifts of electronic surfaces in attached dimers. Such 
shifts are common for alkali atoms and molecules on He droplets |27]. Spectra from 
attached species are shifted relative to what one expects from gas-phase potential energy 
curves [52l ESj [Ml |55l |56] . For certain weakly coupled species, spectra are only shifted 
by a few wavenumbers [27]. Denoting the shift of surface i with Aj, the full state vector 
is propagated according to 



Both surfaces A^S+ and 2^Hg may be affected. Only relative shifts are relevant for the 
signal, such that we cannot differentiate between a shift of the electronic ground - or 
first excited state. Also, a shift of the ionic surface would change the energy distribution 
of the ejected electrons, but not the overall calculated ion yield. This is because we sum 
over all electronic energies in eq. (ig. We find that even large, but fixed shifts of up to 
±100 cm~^ have a very small influence and resemble the gas phase calculation. 

However, it is reasonable to assume that the shifts Aj fluctuate for the following 
reason: The number of He atoms of a droplet is not fixed but varies according to a 




(7) 
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log-normal distribution [5^. Consequently, we have to average over a distribution of 
shifts Aj [58]. Also, for the considered laser intensities, electronic surfaces are slightly 
Stark-shifted. Depending on their position in the laser beam, dimers are exposed to a 
distribution of laser intensities [59]. The beam width leads to a distribution of Stark 
shifts, which in turn has to be treated as random distribution of electronic shifts A,. 

Considering fluctuating potential energies, eq. ^ represents only a single 
realization |\l'j(t)) with shifts Aj(j) in the ensemble. The pump-probe signal is 
proportional to the ensemble average 



where N realizations have been taken into account. In this expression, the final state 
\ipj,3) is obtained from propagation of the realization in which Aj(j) is chosen 

randomly. Through the ensemble average in eq. ([s]), one obtains an incoherent mixture 
of electronic contributions in the final state. The spectrogram of {S{t)) nearly perfectly 
resembles the incoherent sum of contributions, see eq. ^ and fig. |4](c). It does not 
contain any electronic coherence, such that the frequency components Ui are clearly 
resolved. For full decoherence fig. ^c), random shifts in the range of ±5cm~^ around 
the average shifts Ai = Ocm~^ and A2 = — 50cm~^ are sufficient. Note that this 
decoherence is due to the inhomogeneous size distribution and not due to entanglement 
between system and environment [201 EI] • 

To conclude, randomly distributed energy shifts can explain the visibility of the 
ground state component Ux- However, this effect lacks an explanation for the decay 
of the frequency component ua of the excited state, which is observed in the HENDI 
experiment. Therefore, we consider vibrational damping next. 

3.2. Damped vibrational wave packets 

Dissipation can be treated fully quantum mechanically by using approximate master 
equations for the density operator of the (reduced) system. The overall dissipation 
of an excited WP originates in our case from the interaction of vibrational degrees of 
freedom of the molecule with collective degrees of freedom of the helium droplet. Note 
that in recent studies [22] vibrational relaxation rates for alkali dimers on ^He clusters 
are estimated from full quantum Monte Carlo calculations for a few He atoms. These 
rates turn out to be roughly of the same order of magnitude as our phenomenologically 
chosen damping rates below. 

In our approach, the density p of the damped WP in an electronic state \i) evolves 
according to the master equation 



which is of Lindblad form [62]. It describes friction for near-harmonic systems at 

effectively zero temperature in the rotating wave approximation (see, for instance. 




(8) 




(9) 
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Figure 5. (a) Shown is single realization of a wave function in position space. 
The WP circulates in the A state and is damped with a rate 71 — 0.3/ps. (b) 
Coordinate expectation value for that single realization. The noise is due to the 
stochastic propagation, (c) Potential curves A (black) and 211 (red). The latter is 
shifted by one photon energy (A = 800 nm) to clarify possible transitions to the ionic 
state (arrows). A fully damped WP can be mapped to the final state. 



[55]). The first term in eq. (|9| contains the molecular Hamiltonian Hi of a specific 
electronic state \i) and determines the unitary evolution of the WP Relaxation 
of the vibrational WP is achieved through the second, irreversible contribution. Any 
initial state approaches the ground state on a time scale 7"^, i. e. 7 denotes the damping 
rate. In eq. (j9|, a, a''' are the creation/annihilation operators of a harmonic oscillator, 
defined through 

a = ^( J—X + t—^=p] . (10) 
V2 Y V mue ^Jhmuje J 

Here, X and P are the usual position and momentum operator w.r.t. the harmonic 
oscillator minimum and Ue its frequency. The quantum optical master equation ^ 
can be derived from a von-Neumann equation for the full system and is valid only 
for weak couplings between system and "bath" (which is the helium droplet here). 
Also, in the derivation one makes use of the Markov and rotating wave approximation. 
For a significant temperature, additional terms that describe thermal excitations from 
the environment have to be taken into account [51]. The damping constant 7 may in 
principle be derived from a microscopic description of interaction between system and 
bath (Fermi's Golden Rule). We do not specify this interaction, but use the damping 
rate 7 as a fit parameter to obtain agreement with experimental data. The master 
equation ^ induces the evolution of a pure initial state into a state mixture. 

For the numerical propagation, we return from the density matrix description to a 
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Schrodinger-type equation for the state vector. It is not possible to evolve a pure state 
into a mixed state with a deterministic Schrodinger equation. One therefore considers 
a stochastic differential equation for a state vector, quantum state diffusion (QSD) \65\ . 
The density is recovered from the average over several realizations of state vectors. Given 
a master equation in "Lindblad form" , as in our case, it is straightforward to state the 
corresponding quantum state diffusion (QSD) Ito stochastic Schrodinger equation for a 
state vector 



mt)) 



)dt 



[ID 



+ -7 (2(at)a - a^a - |(a)|2) \^)dt + ^{a - {a))\iP) dUf) . 



The left hand side means Idipit)) = lipit + dt)) — \ilj(t)), i. e. the change of the state after 
a time increment dt. The second term induces transitions to lower lying vibrational 
states. The third term is stochastic and contains complex normalized Ito increments 
d$,{t), which satisfy 

de = (dcy = 

d^dC = dt. (12) 

The increment of the density up to second order in dt is given through dp = \dtp){tp\ + 
\'ilj){dilj\ + \dilj){dip\. Since eq. (11) is written in Ito form, one can easily take the average 
with respect to (12) and prove the equivalence with the given master equation ([9]). 
Recovering the density p = \il!{t)){ilj{t)\ from several realizations of {ip) amounts to 
obtaining expectation values through (A) = Ti{Ap) = {ip\A\ip). In fig. sjwe show the 
coordinate expectation value (R), obtained from a single realization of The norm 
of the state in eq. (11) is conserved, i.e. d{{'ip\'ip)) = 0. However, due to the finite time 
step At, the norm can slightly fluctuate. Therefore, for the numerics, we impose norm 
preservation by renormalizing the state vector after every time step. After replacing 
the ladder operators in eq. (11) through their definition in (10), the r. h. s. is strictly 
separable in operators, which act in either momentum or coordinate space. Therefore, 
the split operator method [49J for the propagation of the WP can still be used, which 
is an advantage of this approach. 

Helium induced damping of vibrational wave packets on a single electronic surface 
is described through eq. (11), which is equivalent to the master equation eq. ^ on 
average. In order to obtain damping on several electronic surfaces, we propagate a full 
state vector |^l/(t)) according to 



\d^{t)) 



h 



(13) 



coupling to He bath 

The average is taken over several realizations of state vectors to recover the density via 
p{t) = |\I'(t))(\I'(t)|. In eq. (13), D{jj) is the generalization of the r.h.s. of eq. (11) to 
obtain damping and accompanying fluctuations on an arbitrary surface. 



\lo (2(4)«i - - K«j)l') + - {cij))dij{t) (14) 
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Figure 6. Spectra J^{uj,t) in tlie time-frequency domain at A = 833 nm. (a) Full 
damping, (b) Damping and state-dependent desorption. Also shown is the frequency 
upon averaging in the lined frequency interval. Explanation see text. 



In eq. (14), all stochastic differential Wiener increments dC,j are taken independently 
of each other. In the ladder operator of surface j, one uses the the position w.r.t. 
the respective harmonic oscillator minimum and its frequency. For most parts of the 
following, for simplicity, we set damping constants to be equal jj = 7. However, 



we hasten to add that an interesting exception is provided in section 3.4 The 



agreement with experiment improves, if one allows for undamped vibrational motion 
in the electronic ground state, 70 = 0. 

The damping model provides an explanation for the signal decrease at A = 833 nm 
(excitation scheme I). After several periods, vibrationally damped WPs on the A surface 
no longer enter the initial FC region and therefore the ion yield decreases. In fig. [6] we 
depict the result of the damping model eq. (13), using 71 = 0.15/ps. Clearly visible is a 
shift of the central frequency ua oj'a > ^a, since the WP relaxes to lower vibrational 
levels. A frequency shift is also observed experimentally [1]. In the full damping model, 
however, the signal decays to zero, since at this wavelength vibrational relaxation leads 
to a complete "closing" of the initial FC window. This result is in contrast to the 
experimental observation, where a pronounced oscillation is also present at later delay 
times, see [I] and fig. |2](b), and requires further studies below. 

At A = 800 nm (excitation scheme II), full damping leads to the result shown in 
fig. ^a). Through dissipation, the WP in the state A^S+ leaves the initial FC region. 
However, the decelerated WP approaches another FC window around the equilibrium 
distance, cf. fig. |5} after several circulations. The signal therefore shows a massive 
increase of the (shifted) frequency component u'^. 

For both excitation schemes, the inclusion of damping improves agreement with 
experiment for the first ~ 10 ps. However, at later times, there are still significant 
discrepancies. These can be removed by taking into account desorption of dimers off 
the droplet. 
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Figure 7. Calculated spectra in the time-frequency domain at A = 800 nm. (a) Full 
damping: The WP in the A state leaves the initial transition region, consequently the 
component uja decreases. The fully damped WP, however, reaches a transition region 
at the equilibrium distance, such that uja returns. See also fig. [sj (b) State-dependent 
desorption model, (c) State-dependent desorption model together with undamped 
motion in the ground state. 



The clear visibility of both frequency components u'^ and ux is again attributed 
to electronic decoherence: Since the damping scheme is carried out (stochastically) 
independent on all involved electronic states, damping leads to electronic decoherence 
similar to the fluctuating shifts we assumed previously. The final signal is an incoherent 
mixture of electronic contributions, as previously considered through eq. (|6]). However, a 
massive increase of the component co'^ is not observed experimentally (compare fig. |3](b) 
and fig. |7[a)). 

Desorption of the dimer off the droplet prevents the full vibrational damping in 
the electronic ground or excited state. Indeed, as we will show, desorption implies 
the disappearance of component coa and simultaneous ongoing presence of ux- For 
excitation scheme I, desorption explains the observed signal oscillation at later delay 
times. 

3.3. Desorption 



The release of a dimer into the gas phase takes place at some random time t'. In our 
model, a state vector evolves according to eq. (13) up to time t'. As the He influence 
vanishes at t', shift/damping terms in eq. (13) are set to zero. We have to consider 
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the ion yield at delay time r, which is changed because of the dissapearance of the He 
influence at time t'. If the desorption occurs before the decay of the probe pulse {t' < r), 
the resulting signal is denoted as S'(r, t'). If the desorption occurs after decay of the 
probe pulse (t' > r), the ion yield is unaffected by the desorption process. Anything 
that happens after the probe pulse, will not be mapped to the ion yield. The resulting 
signal is denoted with S{t, t). 

Note that in an ensemble of attached dimers, the desorption time t' will be 
distributed according to some probability distribution P(t'). We therefore have to 
calculate the pump-probe signal for various desorption times t' and consider the averaged 
signal. 

We first assume that the electronic state occupation is not relevant in the desorption 
process. A dimer therefore has a certain probability to stay on or leave (still being in the 
superposition) the droplet, but this probability is independent of the specific electronic 
occupation. We assume a constant probability 

Pos{t',t' + At') = At'RD (15) 

that the dimer leaves the droplet within a small time interval {t',t' + At). Here, Rd 
marks the constant desorption rate. The probability to find the dimer attached at 
arbitrary time t' is hence given through 

Pon(t') = e-*'^"- (16) 
The signal contains weighted contributions and is obtained from 

POO 

{SiT))=RD Ponit')SiT,t')dt' (17) 

^0 

= Pon(r)5(r, t)+Rd [ Pon{t')S{T, t')dt'. 

Jo 

In the second line we have used that for t' > r (desorption after decay of the probe 
pulse) the ion yield does not change. In this case, the dimer is fully damped until the 
probe pulse has passed. 

In a previous study, we made use of the state-independent desorption scheme 



eq. (17), see Ref. [2]. Although it is possible to find reasonable agreement with 
experimental findings, it is an oversimplification to not take into account the electronic 
state for desorption. It is likely that only electronically excited dimers leave the droplet. 
The electronic excitation implies a larger degree of distortion of the helium enviroment 
[H]. On the other hand, it is reasonable to assume that (slowly moving) ground state 
systems stay attached. Therefore, we consider an alternative desorption scheme and 
only allow excited dimers to leave the droplet. Desorption is again described in terms 
of a constant in time desorption rate Rd- The probability for the molecule to desorb, 
however, is now proportional to the excitation probability pe = J2i=/=oPiy such that 



eq. ( 15 ) is replaced by 

Pos{t\t' + At')=p,RDAt'. (18) 
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Note that in this scheme, the state after desorption does not contain any electronic 
ground state component (see later). Those dimers that do not desorb evolve according 



to the dissipative dynamics eq. (13). In order to compensate for apparent loss of 
ground state dimers, we need a third possible channel: With probability PgRo^t' = 
(1 — PejRD^t' the dimer remains on the droplet and is projected onto its electronic 
ground state. For the full averaged signal, we obtain 

(^SD(r)) = Pon{r)S{T, t) + p,Ro I Pon{t')S,{T, t')dt' + (19) 

Jo 



PgRo / Pon{t')Sg{T,t')dt' . 

Jo 



Se{T,t') is the signal obtained upon removal of the helium influence at time t' ^ t and 
subsequent projection on the excited superposition of electronic state. Likewise, 5'g(r, t') 
is obtained upon projection on the ground state. Note that a renormalization of the full 
wavefunction |\E') is required after projection, such that = 1 is always valid. 

In fig. ^h), we show the spectrogram for (^^^(r)) at A = 833 nm (scheme I). 
We use a desorption rate Rd = 0.1/ ps, while the damping/shift parameters are not 
changed (as before , 7 = 0.15/ps and A2 = — 50cm~^). In the spectrogram, a (small) 
frequency shift is ua is observable. The shift is mainly due to contributions 



of vibrationally damped dimers by means of the first term in eq. (19). We obtain 



a frequency shift ua (^'a the model, but the shift is more pronounced in the 



experiment [1]. The second term in eq. (17) marks contributions from dimers which 



are damped up to time t' < r, but are not damped afterwards, i. e. upon release from 
the droplet. Early desorbing dimers are not vibrationally relaxed and the WP in the 
A state continous to reach the initial FC region. From this point of view, undamped 
dimers may well contribute to the ion yield. Indeed, contributions to the ion yield at 
later delay times t > 20ps are exclusively attributed to these undamped dimers. There, 
the oscillation frequency is near the initial gas phase value ua = <^v,v+i- 

As a final note, we find that the observed frequency shift ua oj'a becomes 
negligible if we do not include the electronic shift A2 in the model. 

As the laser frequency is further increased, the vibrational energy of the WP in the 
A^E+ state increases (scheme II). Due to the larger elongation and faster dynamics, one 
may think of larger damping of vibrational motion and/or faster desorption of dimers 
off the droplet. In fig. ^h) we show the spectrogram at A = 800 nm, as obtained 
for 7j = 0.15/ps and Re = 0.5/ps, i. e. we assume a faster desorption as before for 
scheme I. In fact, this value for the desorption rate means that dimers quickly leave 
the droplet after the pulse excitation. Until desorption, electronically excited dimers 
are fully damped. During that time, they approach an intermediate transition region 
with smaller overlap with higher electronic states, see fig. |5|^c). As a consequence, the 
component ua fades away after several picoseconds. Also, vibrational WPs of desorbed 
dimers do not reach the FC window at the equilibrium distance, which is also marked 
in fig. [sj^c). Therefore, an increase of the component ua is excluded, compare fig. [Tj^a) 
and (b). To conclude, damping in connection with fast desorption of dimers explains 
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the experimental result, see fig. [3[b), where the component 00 a is only visible in the 
beginning of the measurement and then disappears. Note that dimers which remain 
on the droplet in the ground state, are fully damped. Therefore, the ground state 
component ux decreases in this model. 

3.4- Undamped ground state WPs 

Agreement with experiment can be further improved, if one allows for undamped 
motion of the vibrational WP in the ground state. This is of particular relevance 
for excitation scheme II, for which the result upon leaving the model parameters for 
damping, desorption and shift unchanged, but setting 70 = 0, is shown in fig. ^c). 
Frequency components at later delay times are exclusively attributed to the ground 
state motion, i. e. the component ux is clearly visible. Undamped vibrational wave 
packet motion has been discussed in terms of a critical Landau velocity i^crit in Ref. [2]. 
The existence of fcrit in the superfiuid nanodroplet may allow for frictionless motion 
of slowly moving ground state WPs. In fig. ^c) the ground state motion is assumed 
frictionless; good agreement with the experimental result, see fig. [sj^b), is also obtained 
for nearly frictionless motion 70 ~ 0. 

Note that for the observation of the ground state WP, the electronic shift A2 is less 
important. It only leads to a slightly different ratio between the ground-and excited 
state frequency component. 

4. Conclusions 

We consider vibrational wave packet dynamics of dimers attached to He nanodroplets. 
It is found that (calculated) gas phase spectra and spectra from dimers attached to He 
droplets are markedly different. The interaction between droplet and dimer influences 
the vibrational dimer dynamics in three ways: Shifts, damping, and desorption. All 
three ingredients are taken into account in a phenomenological manner and each 
contributes to characteristic changes of the resulting spectra. 

First, we study electronic decoherence, occurring for instance due to slightly 
fluctuating shifts of electronic surfaces. We find that, indeed, resulting spectra do not 
show electronic interferences such that contributions from several electronic states, in 
particular the ground state, are clearly resolved. In this way, the modeled spectrum is 
already similar to the experimental finding. 

However, pure electronic decoherence cannot explain a decreasing contribution 
to the signal from wave packets in the first excited A^S+ state. Consequently, 
vibrational damping of wave packets is taken into account, which improves agreement 
with experiment for short delay times. We find that damping is not present over the 
full observation timescale, probably due to desorption of dimers from the droplet. We 
use a desorption scheme, which takes into account the occupation of electronic levels. 
We find that the desorption rate depends on the mean vibrational energy of the wave 
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packet. At A = 800 nm, desorption is fast, taking place on average several picosecond 
after excitation (i?^) = 0.5/ps). 

Note that in this study, all involved electronic states are spin singlet states. For the 
WP in these singlet states, we find a relaxation rate which is significantly higher than in 
previously considered spin triplet Rb2 dimers attached to He droplets. The smaller rate 
in the latter may be ascribed to the orientation of the dialkali axis relative to the droplet 
surface. Recent calculations show that the axis of spin triplet dimers is oriented parallel 
to the droplet surface while singlet states are assumed to be oriented perpendicularly 
[661 [29]. Upon the perpendicular orientation in our case, the dimer might interact with 
the droplet more efficiently, such that relaxation is faster. 

Damping appears to be absent for (slowly moving) wave packets in the electronic 
ground state. There, vibrational motion is found to be nearly frictionless. A potential 
energy surface for the dimer-droplet system, which is currently underway for Rb2 [66] . 
should give additional insight on dissipation rates of the wave packet. 
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